Realizing Dili Class Topological Superconductors using d x 2_ y 2-wave Superconductors 
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In this work, we show that a quasi-one-dimensional 4 2 - y 2_wave superconductor with Rashba 
spin-orbit coupling is a topological superconductor (TS). This time-reversal invariant Dili class TS 
supports two topologically protected zero energy Majorana fermions at each end of the system. In 
contrast to proposals using s-wave superconductors [1HT3] in which a strong magnetic field and the 
fine tuning of the chemical potential are needed to create Majorana fermions, in our proposal, the 
topologically non-trivial regime can be reached in the absence of a magnetic field and in a wide range 
of chemical potential. Experimental signatures and realizations of the proposed superconducting 
state are discussed. 



Introduction — A Majorana fermion is a real fermion 
which has only half the degrees of freedom of a usual 
Dirac fermion. It was first pointed out by Read and 
Green pQ that zero energy Majorana fermion modes ex- 
ist at the vortex cores of a 2D p x + ip y superconductor 
and these Majorana fermions are non-Abelian particles 
[2]. Soon after, Kitaev constructed a spinless fermion 
model and showed that a single Majorana fermion ex- 
ists at each end of a p-wave superconducting wire [3]. 
Recently, several groups proposed that effective p-wave 
superconductors can be realized when an s-wave pairing 
is induced in systems with spin-orbit coupling [4HT3]. 

Particularly, the (quasi)-one-dimensional effective p- 
wave superconductors attracted much attention [6HT3] 
due to the fact that Majorana fermion end states can 
exist in the absence of vortices and the energy separa- 
tion between the Majorana zero energy mode and other 
finite energy fermionic modes is relatively large, on the 
order of the p-wave pairing gap [10] . 

The existence of Majorana fermions in the above men- 
tioned systems is profound. It is related to the symmetry 
of the Hamiltonian which decribes the system. Accord- 
ing to symmetry classification of Hamiltonians [T4HT6] . 
a BdG Hamiltonian with particle-hole symmetry, broken 
time-reversal symmetry and broken SU (2) spin rotation 
symmetry, falls into the D class. In one spatial dimen- 
sion, a D class Hamiltonian is classified by a Z 2 topolog- 
ical number. A system described by a BdG Hamiltonian 
with a non-trivial Z 2 topological number possesses Ma- 
jorana end states. 

To be in the proposed topologically non- trivial regime, 
it requires a Rashba spin-orbit coupling to break the spin 
degeneracy, a magnetic field to break the Krammers de- 
generacy at the Rashba-band crossing point (RCP), fine 
tuning the chemical potential to the RCP and finally 
induce an s-wave superconducting pairing at the Fermi 
energy. However, tuning the chemical potential to the 
RCP, which is near the electronic band bottom, reduces 
the electron density severely and electrons can be easily 
localized by disorder in this regime. The strong mag- 
netic field required can also suppress superconductivity. 
A schematic picture of this proposal is shown in Fig. 2a. 

A Dili class Hamiltonian respects both time-reversal 
and particle-hole symmetry and breaks SU{2) spin- 
rotation symmetry [T4HT6] . In this Letter, we point out 



that a Dili class TS can be realized when electrons in a 
quasi-one-dimensional wire with spin-orbit coupling ac- 
quire a (i a .2_ 2y 2-wave pairing. In the topologically non- 
trivial regime, two zero energy Majorana fermion modes 
appear at each end of the wire. In our proposal, Majo- 
rana fermions can be created in the absence of an external 
magnetic field and in a wide range of chemical potential, 
e.g. there is no need to tune the chemical potential to the 
RCP. In the presence of an external magnetic field, the 
system is in the D class and a single Majorana end state 
appears at each end of the wire. Experimental signatures 
and realizations of this Dili class TS will be discussed at 
the end. 

Single- channel model — Before studying the more re- 
alistic quasi-one-dimensional quantum wires, in this sec- 
tion, we first consider a strictly one-dimensional Dili 
class Hamiltonian which supports double Majorana end 
states in the absence of an external magnetic field. We 
construct the following model 

Hid = H t + Hso + Hsc + H z 

H t = Ej,a -|(^]+l, a ^ia + h.C.) ~ ^)^ja 

H SO = E J>)/3 -|^]+l, a (^)a,/3^,/3 + h.C. (1) 

Hsc = Ej iA (^t +1|t ^ 4 - ^ +1 ,^], t ) + h.c. 

where ipj is a fermion operator at site j, a and f3 are the 
spin indices, t is the hopping amplitude, ctR is the spin- 
orbit coupling strength, Aq is the superconducting pair- 
ing amplitude, and cr y is a Pauli spin matrix. V z denotes 
the strength of the Zeeman term. Without the pairing 
terms, the above model describes a wire with spin-orbit 
coupling. If an s-wave (on-site superconducting) pair- 
ing is induced on the wire, as it is done in Ref . [6HT3] . 
no Majorana fermions can be created without breaking 
time-reversal symmetry. In the following, we show that 
our model supports double Majorana end states in the 
presence of time-reversal symmetry. 

The energy spectrum of H 1 d with V z = is shown 
in Fig. la. Due to Krammers degeneracy, every energy 
level in Fig. la is doubly degenerate. It is evident from 
the energy spectrum that zero energy modes exist when 
the chemical potential satisfies \fi\ < |or|. The sum of 
the two ground-state wavefunctions is shown in Fig. lb 
to comfirm that the zero energy modes are end states. 



FIG. 1: a) Excitation energy versus chemical potential. The 
parameters of Hid are: L = 600, t = 12, Ao = 2 and ctR = 4, 
where L is the number of sites. Zero-energy modes exist in the 
region bounded by the vertical dotted lines where < \(*r\. 
All the bands shown in the figure are doubly degenerate due 
to time-reversal symmetry, b) |^i| 2 + 1^2 1 2 versus x where x 
is the site label. The ground-state wavefunctions are localized 
at the edge in the topologically non-trivial regime \i — olr/2. 
c) In the topologically trivial regime, e.g. \i — 1.5a#, the 
ground-state wavefunctions are predominantly in the bulk. 



In other words, there are two Major ana fermions at each 
end of the wire. In the topologically trivial regime where 
M > |<2r|j the ground state wavefunctions are predomi- 
nantly in the bulk as shown in Fig.lc. 

To understand how the Dili topological superconduct- 
ing state is achieved in our model, we note that a general 
criteria for realizing a one-dimensional Dili class TS is 
to have an odd number of negative pairing amplitude at 
the Fermi points with Fermi momentum between and 
7r [17] . We show that this is indeed the case for H\r>. 

In the momentum space, Hamiltonian Hid can be 
written as 

yt 



H 1D (k) = J2k ^ k [-(tcosk + /i)<7 + a R smka y ]^ k - 
Aocosfc^L^lfci + h.c 



where ^ 



W4t>^jL)- The Hamiltonian has spectrum 



E(k) = ± yj {{—t cos k — fj) ± aR sin k) 2 + (Ao cos k) 2 
and it is generally gapped unless = aR at which 
points topological phase transitions take place. In the 
basis which diagonalize the Rashba term, the Hamilto- 
nian can be written as 

HiD(k) = Efe,a=± [-(tcosk + fi) + a\a R smk\]$l a 'ip ka -\- 
sgn(fc) A cos k$l a ft_ ka + h.c, 

( 3 ) 

where ^) ka denotes a fermion in the new band basis. 
When < |or|, there are two Fermi points fci,&2 with 
< fci, /u2 < 7r. In this regime, it can be shown that one 
and only one of the pairing amplitudes of the two bands 
at the Fermi level Ao cos k\ and Ao cos k<i is always neg- 
ative. Therefore, the superconductor is in the topologi- 
cally non- trivial regime. A schemetic picture is shown in 
Fig. 2b. When > |or|, we have cos&icos/^ > and 
the pairing amplitudes at the two Fermi points have the 
same sign and the system is in the topologically trivial 
regime. 



FIG. 2: a) The energy versus momentum of a typical one- 
dimensional Rashba band with a finite Zeeman term. In s- 
wave-pairing proposals, a strong magnetic field is needed to 
break the Krammers degenercy at the RCP. The chemical 
potential needs to be tuned to the topologically non-trivial 
regime which is bounded by the horizontal dashed lines near 
the electronic band bottom. Electron density in this regime is 
low. b) In our one-dimensional model with V z = 0, topologi- 
cally non-trivial regime is located at the middle of the band 
which is bounded by the dashed lines. In this regime, the 
pairing amplitudes A cos k\ and A cos at the two Fermi 
points, ki and fe, have opposite signs. 



In short, in order to reach the topologically non-trivial 
regime, we need to break the spin degeneracy by the 
Rashba term and induce a ^-dependent pairing such that 
there can be an odd number of negative pairing ampli- 
tudes for positive k at the Fermi energy. It is important 
to note that there is no need to tune the chemical poten- 
tial to the RCP which is near the band bottom. If the 
induced pairing is s-wave [6UT3]. the topologically non- 
trivial phase is not accessible. 

The topologically non-trivial state can be further ver- 
ified by calculating the topological invariant of Hio{k) 
[T7] , The topological invariant can be written as 



( 2 ) N Dili 



Pf[Tq(k = tt)] 
Pf[Tq(k = 0)] 



exp{- 



dkTT[q\k)d k q(k)]} : 



where Pf denotes the Pfaffian 

reversal operator, and q(k) = | [e ze; - ^ (<r — a y ) + 



(4) 

T = icfy is the time- 

1 \„i6-(k)< 



e^+O) (ctq + ct^)] which is an off-diagonal block of the fiat 
band Hamiltonian [T71 [18] derived from Hio(k). Here. 



i6 ± 



-t cos(fe) — /izbcKjj sin(fe)+iAo cos(fe) 



From Eq 



y/[—t cos(/c)— /x±oi.r sin(/c)] 2 + [Ao cos(fc)] 

the topological invariant number can be found to Be 
trivial (Njjjjj = 1) when |/i| > aR and non-trivial 
{Ndiii = -1) when < a R . 

Single- channel model with finite V z — When V z is fi- 
nite, the energy spectrum of Hio(k) becomes E(k) = 

±^F(k)±2^G{k), where F(k) = (tcosk + /i) 2 + 

a 2 R sin 2 k + A 2 , cos 2 k + V 2 and G(k) = (t cos k + n) 2 V 2 + 
(t cos k + ^) 2 c^r sin 2 k + V^ 2 Aq cos 2 k. From the energy 
spectrum, we note that the e nergy ga p closes when 
(fi ± tf = V z 2 - A 2 and | M | = VK 2 + < 

Moreover, the V z term breaks time-reversal symme- 
try and changes the Hamiltonian from Dili class to D 
class. It is known that a ID Hamiltonian in D class 
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FIG. 3: Excitation energy as a function of chemical potential. 
The parameters of Hid are: t — 12, Ao = 1, cur — 4 and 
V z = 10. Points di, bi and a (with i=l,2) denote gap closing 
points with chemical potential \i — =p| yV? — A§ + t\, fi = 

T|\A? + a R\ an( i M = T|\A? - Ag— 1| respectively. Double 
Majorana fermion end states exist when /i is between ci and 
C2. Single Majorana fermion end states exist when fi is in the 
regions between a\ and ci, C2 and a2. 



may support single Majorana end states [3j [14]. The 
energy eigenvalues of Hid versus the chemical potential 
are shown in Fig. 3. In Fig. 3, the region between c\ and 
C2 allows double Maj orana ferm ions where c\ and C2 are 
points with \i = T|\/^ 2 — — ^1 respectively [19]. In 
the regions between a\ and ci, C2 and <22, single Majo- 
rana fermion end state s emerge . Here, a\ and a2 indicate 
points with \i = TlVK 2 — + ^1 respectively. When 
|/i| > \\/V 2 — Aq + t|, the Hamiltonian is topologically 
trivial and the Majorana end state s disappear. At points 
&i and 62 where /i = T\/K 2 + the energy gap closes 
but there are no topological phase transitions at these 
points. 

The even and odd number of Majorana end states in 
Fig. 3 can be verified by calculating the Z2 Majorana 
number M of Hi D (k). Following Refs [3] [11] the Ma- 
jorana fermion number can be defined as 



M = sgn[Pf£(0)]sgn[Pf£(7r)] = ±1. 



(5) 



The matrix B(k) is defined as B(k) = HiD(k)(a x (g) ao). 
B(k) is anti-symmetric and its PfafHan is well-defined. 
A4 = ±1 indicates the even and odd number of Majo- 
rana fermions at one end of the wire respectively. In 
terms of the parameters of the Hamiltonian, the Majo- 
rana number can be written as 

sgn{ [(* + - (V? - Aq)] [(-* + n) 2 - (V? - Aq)]}. (6) 

The Majorana numbers calculated according to Eq(6]are 
consistent with the results shown in Fig. 3. The Majo- 
rana number is A4 = — 1 when \\/V 2 — Aq — t\ < |/i| < 
I y/V? - A 2 , + 1\ and M = 1 otherwise. 

Multi-channel case — In realistic situations, multiple 
transverse sub-bands of a wire are occupied and it is im- 
portant to show that Majorana fermions exist in this sit- 
uation. In this section, we show that Majorana fermions 
exist in quasi-one-dimensional wires. Importantly, the 
quasi-one-dimensional model can be realized by induc- 
ing d x 2_ y 2-wdive superconductivity on a wire with strong 
spin-orbit coupling. 




Majo rana fermjons 
Rashba system 
uperconductor 



FIG. 4: A wire with strong spin-orbit coupling in proxim- 
ity to a d cc 2_ y 2-wave superconductor. Double Majorana end 
states may appear in the absence of an external magnetic 
field. 



In the quasi-one-dimensional case, the Hamiltonian 
can be written as: 

H q iD = H t + Hso + Hsc + Hz, 

H t = E R ,d,c* -|*(^R+d,a^Ra + h ' C ') ~ ^R,a^Ra 
H SO = E R ,d,c^ "I a ^R +d , a Z • (°aP X d)^R,/3 + ft.C. 
H SC = E R |[ A 0(^R +dx ,t^R,; - ^R+d x ,^^R,t)~ 

A o(^+d t^kl ~ ^R+d y ,^R,t) + h ' C '} 
H Z = Er Vz (^R t ^Rt - ^r;^R|) • 

(7) 

Here, R denotes the lattice sites, d denotes the two 
unit vectors d x and d y which connects the nearest neigh- 
bor sites in the x and y directions respectively. This 
model is the same as the tight-binding model in Ref[T2] 
except for the superconducting pairing terms. The pair- 
ing terms in H q \D can be written as A [cos(/c :c ) — cos(fc y )] 
in the momentum space. Therefore, H q \D describes a 
quantum wire with spin-orbit coupling and a d x i _ y 2 -wave 
superconducting pairing. A schematic picture of the ex- 
perimental setup is shown in Fig. 4. 

The energy spectrum of the Hamiltonian with V z = 
is shown in Fig. 5a. The length of the wire is chosen 
to be much larger than the superconducting coherence 
length L ^> t/Ao and the width is comparable to the 
coherence length W « t/A$. In sharp contrast to the 
s- wave-pairing proposals in which the topologically non- 
trivial regime can be reached only when the chemical 
potential is near the RCP, in our proposal, zero energy 
Majorana modes live all over the full band even when 
the Fermi level is far away from the band bottom. It is 
also evident from Fig. 5a that the zero energy Majorana 
modes are separated from other fermionic states by a 
large energy gap which is of order Ao. Due to time- 
reversal symmetry, all energy levels in Fig. 5a is doubly 
degenerate. Therefore, in the topologically non-trivial 
regime, there are two Majorana fermions at each end of 
the wire. It can be shown that these double Majorana 
end states are robust to disorder. 

In the presence of the V z term, time-reversal symmetry 
is broken and the Hamiltonian is in D class. To break 
some incidental deg eneracy, a ra ndom on-site potential 
V with strength ^ (V 2 ) - (V) 2 = A /2 is introduced 
into the system. The resulting energy spectrum of the 
model with V z = 3 is shown in Fig. 5b. Without time- 
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FIG. 5: a) Excitation energy versus chemical potential. The 
parameters in H q iD are V z = 0, t = 18, Ao — 2 and aR — 
4. The width of the wire is W — 6 and the length is L — 
600. All the states shown are doubly degenerate. Zero energy 
Majorana modes, denoted by the red lines, can appear in a 
wide range of chemical potential. The normal state electronic 
band bottom is denoted as BB b) Excitation energy versus 
chemical potential in the presence of random potential. The 
parameters in H q iD is the same as in a) except V z = 3. The 
topologically non-trivial regimes with an odd number of zero 
energy modes are bounded by dotted lines. 

reversal symmetry, the spectrum is no longer doubly de- 
generate. From Fig. 5b, it is evident that zero energy 
modes which represent single Majorana end states ex- 
ist. As expected, the single Majorana modes are more 
robust when the chemical potential is near the band bot- 
tom. Majorana fermions exist when the chemical poten- 
ital is further away from the band bottom, however, the 
mini-gap, which is the energy gap between the Majorana 
modes and other finite energy fermionic modes, is much 
smaller as in the s-wave pairing case [20] . 

It is shown previously that a single Majorana fermion 

2 

induces a zero bias conductance peak of G = 2^- in An- 
dreev reflection experiments [21, 22 when a normal metal 
lead couples to a Majorana end state. On the other hand, 
the double Majorana end states in Dili class TS induce a 
zero bias conductance peak of G — 4^- [23, 24 instead. 
Therefore, the two topological superconducting phases 
can be distinguished from each other and from trivial su- 
perconducting states through Andreev reflection experi- 
ments. 

Discussion — A few important comments follow. First, 
the Dili class topological superconductor proposed in 
this work is truly different from the time-reversal in- 
variant superconducting state obtained by inducing an 
s-wave or a d-wave pairing on the surface of a topological 
insulator (TI) [25H27] . The TI surface state is described 
by a single Dirac cone (in the simplest case). With only 
a single species of fermion at the Fermi level, one cannot 



create a superconducting wire with double Majorana end 
states even though single Majorana fermions can be cre- 
ated in the presence of an external magnetic field. Sato et 
al. studied a d-wave superconductor with Rashba terms 
in the presence of an external magnetic field but no Dili 
class TS phase is reported [28] . 

Second, it is pointed out recently that single Majorana 
fermion end states can be created at arbitary chemical 
potential if a magnetic field is applied along the wire in 
the s-wave pairing case [20j[27]. However, the magnetic 
field B required to reach the topologically non-trivial 
state is still strong comparing to the induced pairing gap 
As, namely, fio\B\ > As where fio is the effective mag- 
netic moment of electrons. 

Third, for simplicity, we assumed that the wire is 
aligned along the x-direction. Tilting the wire with re- 
spect to the x-axis slightly has no major effect on the 
Majorana end states since this kind of perturbation does 
not change the symmtry class of the system. 

Fourth, one possible realization of H q \r> is by induc- 
ing gL^ -wave superconducting pairing on the Au (111) 
surface state. It has been shown that the Au (111) sur- 
face has a Rashba band with Rashba energy of about 
60meV [29]. The proximity gap induced on Au by LSCO 
can reach lOmeV [30]. These large energy scales make Au 
on LSCO a promising candidate of realizing topological 
superconducting states. However, the induced proxim- 
ity pairing shown in the recent experiment may not be 
d-wave due to the presence of strong disorder in the Au 
layer. 

Another candidate material of a Dili class TS is a 
layered heavy fermion superconductor CeCoIn 5 . Bulk 
CeCoIns is a d x 2_ y 2-wdive superconductor. Recently, su- 
perconducting thin films of CeCoIns with only several 
atomic layers thick can be fabricated [51] . Suppose we 
have a thin film of CeCoIn 5 , due to the strong spin-orbit 
coupling and the broken of inversion symmetry on the 
surface, the top layer of CeCoIn 5 acquires a Rashba term 
[32] . Therefore, the top layer of CeCoIn 5 can be de- 
scribed by H q i]j. In a separate work, we show in detail 
that a thin film of CeCoIns is a Dili class TS. 

Conclusion — We show that a quasi-one-dimensional 
<i :c 2_ 2y 2-wave superconductor with Rashba spin-orbit cou- 
pling terms is a Dili class TS which supports double Ma- 
jorana end states. In the presence of a magnetic field, sin- 
gle Majorana end states appear. These two topological 
superconducting states can be probed using Andreev re- 
flection experiments. We suggest that Au wires on LSCO 
and thin films of CeCoIns are candidate materials of this 
new topological superconducting state. 
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